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Abstract 

We analyze how corrections hnear in the effective range, tq, affect quantities in the three-body 
sector within an effective field theory for short-range interactions. We demonstrate that observables 
can be obtained straightforwardly using a perturbative expansion in powers of tq. In particular, we 
show that two linear-in-ro counterterms are needed for renormalization at this order if scattering- 
length-dependent observables are considered. We exemplify the implications of this result using 
various three-body observables. Analytic results for the running of the next-to-leading-order por- 
tion of the three-body force in this effective field theory are provided. Expressions which incorporate 
O(ro) corrections and relate the positions of features observed in three- atom recombination near a 
Feshbach resonance are presented. 
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I. INTRODUCTION 



Symmetries are one of the most important concepts in modern-day physics. They are the 
foundation of the standard model and are also frequently used as a basis for building new 
physical theories. Even if a symmetry is only approximately fulfilled it can, nonetheless, serve 
as a starting point for a systematic description of physical observables, since the symmetry 
breaking effects may be accounted for perturbatively. This approach has been successfully 
employed in various effective field theories (EFTs) which are low-energy expansions in the 
ratio of a small parameter over a large parameter. The small parameter is frequently a 
small momentum and/or energy associated with explicit symmetry breaking. A prominent 
example of such an EFT is the chiral EFT, whose starting point is the chiral limit of QCD. In 
this limit, pions are the Nambu-Goldstone bosons associated with the spontaneously broken 
chiral symmetry of QCD. In the chiral EFT the effect of the nonzero (up, down, and strange) 
quark masses is then systematically included order-by-order in the low-energy expansion. 

While the use of symmetries to constrain theories is usually associated with particle 
physics, or possibly many-body physics, it has also become important in other fields, such 
as non-relativistic few-body physics. One example of an important symmetry in this area is 
discrete scale invariance. Vitaly Efimov showed in 1970 that the non-relativistic three-body 
system displays discrete scale invariance if the two-body scattering length is large and the 
range of the interaction is zero PP . One well-known consequence is that in the limit of infinite 
scattering length the ratio of the binding energies of two successive three-body bound states 
is approximately 515. 

Efimov also pointed out that his results apply to systems where the two-body scattering 
length, a, obeys \a\ ^ i, with i the natural length-scale of the two-body potential [1]. They 
are therefore relevant for a number of systems. One example is nucleon-nucleon scattering, 
where the scattering length is large compared to the range of the internucleon interaction. 
A second is provided by manipulation of the atom-atom scattering length of trapped atoms 
by an external magnetic field such that the atomic scattering displays a Feshbach resonance. 
Few-body systems of nucleons and atoms, all with interactions tuned such that a shallow two- 
body bound state is present, will therefore display discrete scale invariance or the remainders 
of this symmetry. For a recent review of this topic see, e.g., Ref. [2]. 

Efimov's results have now been rederived as the leading-order (LO) prediction of an 
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EFT containing only short-range interactions [Sj H]. At LO this EFT corresponds to the 
Efimovian, large-a, i = hmit. We therefore refer to it hereafter as short-range EFT 
(SREFT). The contribution of each subsequent order in the SREFT expansion is suppressed 
by an additional power of i/a, with the first corrections that have to be included in this 
EFT being those resulting from a finite two-body effective range, tq. The effects of these 
corrections have been analyzed over the last years in a number of works |5H7]. In these 
studies effective-range corrections were considered for systems in which the scattering length 
remains fixed. A full analysis requires, however, that we allow for a variable scattering 
length. Such an analysis was reported in |S1 E]. In this paper we lay out the derivation 
of the results reported in [9]. We discuss the renormalization of the SREFT in the three- 
body sector at next-to-leading order (NLO) in the £/a expansion. We show that if, and 
only if, scattering-length-dependent observables are considered, an additional three-body 
counterterm is required for renormalization at next-to-leading order. The analysis of Ref. [8] 
overlooked this result since the combination of the non-perturbative treatment of effective- 
range corrections and the cutoffs A ^ l/rg employed there modifies the ultraviolet properties 
of the theory. The results of Ref. [8] therefore are strict EFT predictions only in the limit 

In Sec. [TT] we will introduce SREFT and briefly review what is known about the two- 



body sector up to next-to-leading order. We then discuss, in Sec. Ill, the so-called modi- 
fied Skorniakov-Ter-Martirosian integral equation which constitutes the application of the 



SREFT to the three-body sector at leading order. Section IV describes how NLO correc- 
tions affect a variety of three-body observables. Section |V] discusses the running of the 
NLO parts of the three-body force, and presents analytic results for these counterterms' 
renormalization-group evolution. Finally, Sec. |VI| presents relations between observables 
measured in three-atom recombination which incorporate effects that are NLO in SREFT. 
We end with a summary and outlook. 



II. THE EFT FOR SHORT-RANGE INTERACTIONS 

We employ the SREFT that describes non-relativistic particles interacting through a 
finite-range interaction with a large scattering length. The inverse of the range of the inter- 
action sets the breakdown scale for this EFT, which is constructed from contact interactions 
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alone. In nuclear physics this EFT is also known as the pionless EFT (see Refs. [TOt [TT] for 
reviews). However, the SREFT can also describe other systems such as atoms close to a 
Feshbach resonance. Since we perform our analysis for identical bosons we will refer to the 
interacting particles as atoms. 

At the heart of any EFT lies the Lagrangian. It includes all possible interaction terms 
that are allowed by the symmetries of the underlying interaction and is built from fields 
that correspond to the degrees of freedom included in the EFT. The SREFT in its original 
form is therefore built from atom fields alone, however, it has proven useful to employ a 
Lagrangian that contains a dimer field. The relation between these forms of the Lagrangian 
is elucidated in Refs. [H [T2]4T^ . In these works the dimer field is sometimes chosen to be 
static, sometimes dynamic. The latter version is particularly convenient for the inclusion of 
effective range corrections, and so we will work with this form of the SREFT Lagrangian, 
written as: 

C = il)^ (^d^ + T/' + aTt l^z^o + ^ - T- [T^^j + h.c) +hT^Tij^ + .. . , (1) 

where the ellipses represent additional higher-order interactions which are suppressed at low 
momenta. As pointed out in Ref. a positive effective range can only be described by 
this theory if a = — 1. In this case the free theory describes an atom and a ghost field. 
We will employ the Lagrangian ([T]) and will expand all elements of all Feynman graphs in 
powers of tq. The SREFT Lagrangian is then equivalent, order-by-order in the ro expansion, 
to the SREFT Lagrangian with an atom field alone, as presented in, e.g., Refs. [TOl [12]. 

The Feynman rules are derived from Eq. ([T]) and the atom propagator in momentum 
space is 

iS{po,p) = ^, (2) 

where po is the energy and p = |p|. The large scattering length leads to large loop effects, 
and the EFT power counting requires therefore that the two-body interaction is iterated to 
all orders (Fig|T]). The resulting dressed dimer propagator is 

^'^ipo,p) = — I ,7 — • 

p,-l--A + '^^-mpo + ^-^-ze + ze 

If the scattering length a is positive, the two-body system will display a bound state {shallow 
dimer). Equation ([s]) therefore has to have a pole at the on-shell four-momentum pq = 



4m m ' 



- — — where ——is the binding energy of the shallow dimer. We can rewrite this condition 



and obtain 



which has the solution 



2 2 

- 7^ - mA + ^7 = 0, (4) 



mV / / 647r2A\ 

We can relate the coupling constants g and A to scattering length a and effective range rg 
by using the S-wave effective range expansion 

1 1 2 

- - = -7 + -ro7 (6) 

Qj Zi 

which leads to a = and ro = ^^^2- Note that o = —\ was used in Eq. ([T]) to derive 

these results. In terms of these quantities, the wave-function renormalization factor Z for 
this bound state is: 

l = z^(,p(p))-i| , ,2=^(l-ro7). (7) 
Thus, the dressed dimer propagator can be reexpressed as 

iP(po,p) = -lATijmg ^ 

—7 + 2^q{i^ + "^Po — + y'—mpo + p^/A — ie + ie 
Two-body scattering phaseshifts for relative momentum k can be obtained from the 
propagator V{po,p) by evaluating it at the on-shell point of two atoms scattering with 
center-of-mass momentum p = and E = k'^/m, and multiplying it by —ig"^. This yields a 
two-body amplitude: 

An 1 

^ (9) 

m 7 — ^ro(72 -I- k"^) + ik' 

i.e. one in conformity with the effective-range expansion around the two-body bound-state 
pole. 

However, the propagator in Eq. ^ cannot be directly employed in an integral equation 
for three-body physics, since it contains a spurious pole. This deep bound state (pole at 
momentum ~ ^/tq) arises since we have chosen to work with a ghost field that leads to 
unphysical features in the theory, such as negative- norm states and/or violations of unitary. 
However, when we expand the propagator in powers of 'jro, the deep pole does not appear 
at any finite order: 

iV(p„p) = j: (r2)-(,+ ^_,„,„+,V4)". (10) 



FIG. 1: The series that is summed to obtain the dressed dimer propagator, leading to Eq. Q. 
The leading-order (LO) dimer propagator is therefore given by 

^Vi«>(p„p) = -■"^'"'^l . . (11) 

-7 + v -mpo + P"^ - te + te 

and the next-to- leading order O{'yro) correction is 

tV- >{po,p) = -t ^x— =. (12) 



III. THREE-BODY AMPLITUDES 

The calculation of three-body observables requires the solution of an integral equation 
which we will call the modified Skorniakov-Ter-Martirosian equation. It is the Faddeev 
equation for non-relativistic particles interacting through two-body and three-body zero- 
range interactions. We will show this equation below. Details of its derivation can be found 
in Refs. [31 II]- A three-body interaction is included there, to ensure cutoff-independent 
results. The running of the associated coupling constant (h) displays a limit cycle. 

Our analysis concerns the corrections to the modified STM equation and is based on 
rewriting all involved quantities in the form 

V{p) = V^'^\p) + V^^\p) + . . . 
t{q,p) = t(°)(g,p)+t«(9,P) + --- 

H{A) = /fo(A) + ^i(7,A) + ..., (13) 



where H{A) = A h/2mg , and each quantity in Eq. (13) is expanded in powers of kro and 
'jrQ. Hence all quantities in the calculations that follow are computed up to next-to-leading 
order in the i/a expansion. Since 7, ~ 1/a and we assume tq ~ i, this implies that our 
goal in this work is to compute t*-^^ and associated NLO parts of physical observables in this 
expansion. The t-matrix obtained in this way is then related to the renormalized amplitude 
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FIG. 2: The diagrams for the next-to-leading order t-matrix. Propagators and vertices with gray 
squares denote NLO corrections. 



for s-wave atom-dimer scattering at relative momentum k, T{k), via: 

T{k) = y/Z{to{k, k; E) + ti(A;, k-E) + .. .)v^, 



(14) 



where, for elastic atom-dimer scattering, the magnitude of the relative incoming momentum 
|k| equals the relative outgoing momentum |p| and the three-body energy, is given by: 

4m m 

In Eq. (14), Z is the wave-function renormalization factor given in Eq. ([T]). We also expand 
this Z-factor in powers of 7ro: 



(16) 



and so on, thereby generating an expansion for T{k) analogous to those listed in Eq. (13). 

The leading-order three-body amplitude, to, is computed by iterating the two- and three 
body interactions of Eq. ([T]) to all orders via the integral equation |3] 



fo(g,p; E) = M{q,p- E) + 



TC 



dq 



J2 



-7+ v/3g'V4-m^- 



-.M{q',p;E)to{q,q';E), (17) 



where 



1 f q^+p^ + qp-mE \ 2H^{A) 

M{q,p; E) = — log + — 

qp \q'^ + p'^ — qp — mE J 



(18) 



Note that we have rescaled the t-matrix in Eq. (17) tQ(q,p;E) = mg to(q,p; E) such that 



it depends only on physical quantities. The integral equation without the three-body force 
was first derived by Skorniakov and Ter-Martirosian and corresponds to three particles 
interacting through zero-range two-body interactions [22] . 
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The LO three-body force, Hf){A) is fixed by fitting it to one three-body observable. A 
combination of analytical arguments and numerical studies then show that Hq{A) varies 
with A as: 

sin(so ln(A/A) + arctan(so)) 



i^o(A) = c- 



(19) 



sin(so ln(A/A) — arctan(so)) 
The constant c is 0{1) and depends on details of the regularization of the modified STM 



equation. For the regulator employed in Eq. (17) we confirm the value c = 0.879 determined 
numerically in Ref. [21]. 

We calculate corrections to the amplitude obtained from the LO modified STM equation, 
to, by considering diagrams with a single insertion of the NLO dimer propagator. It was 
shown by Hammer and Mehen |5] that such a perturbative inclusion of effective range cor- 
rections also requires the insertion of a subleading, energy- independent, three-body force. In 
Fig. |2] we display the diagrams that have to be evaluated. The application of the Feynman 
rules gives, for the first-order correction to the amplitude 



X^tW(q, q', q',-E + ^) ztW(q', p,E-f- 

2m 2m 

,2; 



+ i 



2mg^Hi{^,A) 



A2 



1 + 



q' 



^ .I)(°)(g^, q') iS{E - g^, q') ^to(q, q', ^lo'E^ — 



1 + / ^'D^'\q'o^ q') ^S{E - q'„ q') zto(q', V,E - ^ 



(20) 



The complete NLO correction to the S-wave projection of the t-matrix is therefore 



ti(g,p; E) 



TT 



dq q 



' '2 



-7 



7 + ^3g'74-mE 
^J'iq'^/A-mE -ie 



to{q,q;E)to{q',p; E) 



2^1 (7, A) 
A2 



2 




dg'- 




f 


TT 


Jo 






2 




n 




f 




TT 


Jo 



j2 



-7 + A/3g'V4 - mE - ie 



to{q,q';E) 



dq - 



j2 



-7 + ^J'iq''^/A - mE - ie 



toiq,p; E) 



,(21) 



where ti and Hi are defined via ti(g,p; E) = rQmg'^ti{q,p] E), and Hi = r^Hi. 

The expression above is similar to the one obtained by Hammer and Mehen in Ref. [5]. 
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However, in their work one contribution from the contour integration in Eq. (20) was erro- 



neously omitted. In that approximation the first term in Eq. (21) becomes: 

I fdq'q'' , =to{q,q';E%{q',p;E). (22) 

TT Jo -7 + A/3g'2/4 - mE - le 

As will become clear below, this omission leads to the conclusion that only one counterterm 
is required to renormalize the three-body problem in SREFT at NLO. While this works 
in practice for fixed-a observables, such as the neutron-deuteron phase shifts computed in 
Ref. it produces incorrect results if experiments in which a varies, like those examining 
loss features in gases of cold atoms, are analyzed. 

to, ti, etc. are generally complex when scattering states are considered, albeit real when 
the bound-state problem is considered. It is therefore convenient to introduce the real K- 
matrix, which contains the same information as the t-matrix but is easier to calculate. The 
LO half-on-shell K-matrix obeys an STM equation in which the ie prescription is replaced 
by a principal-value integration (indicated by V) 

K,{k,p-E)=M{k,p-E) + -V j dq' ^ '-Miq',p;E)Ko{k,q';E). 

(23) 

The half-on-shell t-matrix at LO is related to the K-matrix via 



Uk,p;E)^ ^ (24) 



Mk,p-E) 
/-Ko{k,k;E) 

The fully-off-shell t-matrix at leading order is also related to the K-matrix through a similar 
transformation 

fo(g,p; E) = Ko{q,p; E) + i^Ko{k,p- E%{q, k; E). (25) 
This allows us to write the fully-off-shell K-matrix as 



Ko{q,p;E) = M{q,p;E) + —V dq M{q ,p; E)Ko{q,q ; E), 

671 Jq q — k 

(26) 

with E and k related as in Eq. (15). The half-on-shell NLO K-matrix can then be expressed 
in terms of the half-on-shell LO and NLO t-matrix and the on-shell LO t-matrix, LO K- 
Matrix and NLO K-matrix 

Mk,p; E) = ^.ilff] - t^toiKp; E) (fc, k- E) + ^Ko{k, k; E)] , (27) 
l + i^to{k,k;E) 3 L J 

9 



or vice versa 



t\{k,p; E) 



Ki{k,p-E) 



t'-fKoik,p;E) K,ik,k;E)+^K,ik,k;E) 



l-i^Ko{k,k;E) 



+ 



l-i^Koik,k;E) 



(28) 



At next-to-leading order, the half-on-shell K-matrix is then given by the following principal- 
value integral: 



K,ik,p;E) = -VI ,,V. ^+v/3g'V4--i^ 



7r 



2^1 (7, A) 

A2 



-7 + ^3g'74-mE 

2 

1 + -P dq' 



Ko{k,q';E)Ko{p,q';E) 



-7+ ^3g'2/4 -mE 



Koik,q';E) 



2 

1 + -V dq'- 



TC 



j2 



-7+ v^3g'V4-m£; 



(29) 



IV. THREE-BODY OBSERVABLES AT NLO 

In this section we will discuss how different observables are calculated at NLO in our 
perturbative approach. We will consider not only obvious observables, such as phaseshifts 
and three-body binding energies, but also observables typically measured in experiments 
with ultracold atoms. We follow the strategy outlined above and calculate all quantities as a 
series in powers of 7ro and/or k tq. The leading order in this series is then the universal result, 
and the NLO pieces we will derive here encode the first corrections "beyond universality" . 

A. Phaseshifts 



The amplitude T{k) for atom-dimer scattering is related to the atom-dimer S-wave phase- 
shift through 



T{k) 



m k cot 6{k) — ik 



(30) 



The scattering amplitude T{k) in Eq. (30) can also be expanded in powers of 7ro 

37r 1 



nk) 



m k cot 5q + To [k cot 5] ] 
37r 
m 



1 



I ik 

[k cot 5]i 
kcot6o — ik ^ {k cot 5o — iky 
To{k)+T^{k) + ..., 



+ ... 



(31) 
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where the dots refer to corrections beyond NLO. Here kcotS is expanded as 



k cot 6 = k cot 60 + tq [k cot 5] 1 + . 



(32) 



where the [] 1 indicates the part of k cot 6 that is the coefficient of the order tq term in 
the expansion in powers of tq. At leading order we recover the famihar relation between 
phaseshifts and K-matrix 

3 1 

k cot So = —Ko^\k,k;E), (33) 
87 



but our expansion also leads to a relation for the NLO part, that stems from Eqs. (31) and 



(28): 



[k cot <5]i = -—Ko'\k, k;E){^ + K^{k, k; E)/Ko{k, k;E)). 



(34) 



As long as k cot 6q is not large this will yield a correction of relative size ^vq to the leading 
part of k cot 6. 

B. Bound States 

The t-matrix is real for energies below the scattering threshold. The NLO correction to 
the leading order t-matrix is then given by 



ti{q,p; E) 



1 /"^ , , ,2 7 + V3g'V4 - mE ~ 



TT 



dq q , 

-7 + V3g'74 - mE 



to{q,q']E)to{p,q']E) 



2^1(7, A) 
A2 



1 + 



TT 



A 



dq- 



+ v/3g'V4 - mE 



-7 



TT 



dq - 



(I 



12 



to{q,q';E) 



-to{p, q';E) 



(35) 



10 -7 + V3g'V4 - mE 

where we have explicitly dropped the ie prescription. The full t-matrix has a pole aX E = B 
if a three-body bound state exists with this energy. At leading order this implies: 



^o(g,p;i?) = |^ + 7^o(g,p;i?) 



(36) 



E — Bq 

where the function IZq is a regular part. The residue Zq depends on the incoming and 
outgoing momenta q and p, but not on the three-body energy E. Although, in general, more 



than one bound state exists for the systems under consideration here, the decomposition (36) 



is the appropriate one if we are focusing on the NLO shift for a particular bound state. The 
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fact that LO bound-state energies are separated by a factor as large as 515 allows us to 



employ the decomposition (36) for these purposes 



When considering the NLO correction, we must account for both the pole's position (i.e. 
the three-body binding energy) and the residue being shifted by an amount proportional to 
tq. Therefore we have 



to + roh 



E — Bq — Bi 



+ 7^o + 7^l 



Zn ZnB\ Z\ „ „ 

+ V; so + ^ +7^o + 7^l, 



(37) 



E-B^ (E- 5o)2 E-Bo 
where, as usual, the subscript 1 indicates the parts which are first order in tq. In particular. 



Eqs. (37) and (36) imply that the first-order part of t, ro^i, has a pole of order two at 



E = Bn 



roti{q,p; E) = — — - + — + 7^l(g,p; E). 



(38) 



{E - Boy E-Bo 

The residue of this double pole is then related to the shift in the three-body binding energy 
that is linear in the effective range: 

MmE^BoiE - Bo)Hi{q,p; E) 



Bi = ro- 



Zo{q,p) 



(39) 



Equation (39) seems to indicate that the incoming and outgoing momenta, q and p, affect 



the three-body binding energy shift Bi. However, we would expect that the binding energy 
is independent of the incoming and outgoing momenta. This apparent contradiction can be 
avoided if Zo(g,p) is separable with respect to q and p. Therefore, the residue function is 
defined as 

Zoiq,p) = T{q)Tip). (40) 



By substituting Eqs. (40) and (36) into Eq. (17) and taking the residue at E = Bq we find 

"A „/2 



T{q) = - I dq' Miq,q';Ey 

IT 



=r(g')- 



(41) 



'0 -7 + V3g'V4 - mE 

The function r(g) is thus a solution to a homogeneous integral equation, and the overall 
normalization is not immediately determined. We fix this normalization by the condition: 



r2(g)= hni iE~Bo)ioiq,q;E). 



(42) 
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With r(g) in hand, we can insert Eq. (40) into Eq. (35), and multiply by {E — BqY and 
take the limit as E ^ Bq. This yields a result for Bi that is independent of q and p: 

dq q^ 
8ifi(A)ro 



vr 



7 + V3g74 - mB, 



A 



uiBq 

^2 



-\ 2 







-7 + v/3gV4 - m5| 



:r(g) 



(43) 



And indeed, in numerical calculations, Eqs. (39) and (43) prove to be equivalent. Bi can be 



obtained from Eq. (39) if desired, and the result found in that way is independent of q and 
p. 



C. Three-Body Recombination 

Three-body recombination is a collision process in which three free atoms combine into 
a dimer and an atom. The atoms can either recombine into deeply bound two-body states 
{deep dimers) whose properties cannot be described by the SREFT (but see Ref. [I3] for a 
discussion on how the effects of these deep states can be included in the theory), or, provided 
the scattering length a is positive, into the two-body bound state which is explicitly included 
in SREFT. The energy that is released when the two-atom bound state forms is converted 
into kinetic energy and atom and diatom are lost from the trap. The loss rate of atoms in 
a cold atomic gas due to three-body recombination into the shallow dimer is determined by 
the scattering amplitude for the reaction A + A + A^A + D, as we shall now show. The 
scattering-length dependence of the loss rate therefore provides an experimental signature 
of Efimov physics in trapped systems of ultracold atoms. 

The atom loss rate is expressed as 

dn , ^ . 

where n is the number density of free atoms. (The factor of 3 arises due to the loss of three 
atoms in each recombination event.) According to Fermi's golden rule, 

iy,, = 27r|T(p,)r^, (45) 

where T is the amplitude for three- atom recombination: A + A + A — > A + and the 
density of atom-dimer states dvj is 
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FIG. 3: Three-body recombination at leading order 



The kinetic energy at a momentum pf in the atom-dimer system is 

3p/ 



and so the transition rate becomes 



W 



2m 4m 



4m 



(47) 



fi 



2m 
3^ 



P/|trec(p/)P 



(48) 



The recombination rate a is conventionally defined as 



and so 



a 



dn 
dt 

m 
~ 97r 



—San 



P/l^rec(p/)r. 



(49) 



(50) 



At zero temperature three-body recombination takes place at the three-atom threshold and 
the value of the relative momentum pj in the atom-dimer system is therefore 2'j/\/3. 

The y4 + y4 + y4 — > A + D amplitude at leading order in Fig. [3] is related to the half-on-shell 
atom-dimer scattering t-matrix by 



48712 - / 27 



m 



(51) 



which, combined with Eq. (50), determines the leading-order recombination rate as: 

2 



"0 



5127r' 



m 



to|0,^;0 



(52) 
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FIG. 4: Three-body recombination at next-to-leading order 



Upon introducing the effective-range correction, we obtain the sum of the leading-order 
and next-to- leading-order A + A + A^A + D amplitude (Figji]). This is expressed by 

tiel{Pf) + till{Pf) = 3-{-tV2g)^/Zo + z,- [^p(°)(o,o)to + ^^^^'Ho,o)to + ^^^^°Ho,o)^l] 



487r3/2 



2VZo 



m7 



487r3/2 



m7 



1/2 



7ro4(0, 27/ V3; 0) + rofi(0, 27/^8; 0) 



Therefore, the recombination rate is given at NLO by: 



5127r2 



V3 



m 



fo(0, ^; 0) + 7rot~o(0, 0) + rofi(0, ^; 0) 



+ 



(53) 



(54) 



1. Recombination minimum 



At leading order the recombination rate a is related to Kq by 



~, 2, 



512vr2 i^o70,^;0) 



"0 



m 



1 _ ,-167!. ?r r 27 27 . r|\ 



V3' v^' 



(55) 



Therefore, the recombination minimum is determined by the condition Kq(0, ^;0) = 0, 



which leads to to(0, 0) = 0. 

Because ao = at the minimum, the NLO recombination rate in Eq. (54) becomes Olr^) 
in the vicinity of the leading-order minimum. This means that in what follows we cannot 
neglect terms ~ Tq, which complete the square and guarantee that a > 0. 

In order to calculate the NLO correction to the recombination minimum, we evaluate 



«! from Eq. (54) at 7 = 70 + A70, and expand it in powers of tq. First, we expand the 
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scattering amplitude before we square it: 



(1 + 7ro)to(0, ^; 0) + rofi(0, ^; 0) 



d ~ 27 



7=70 



A7 + rofi(0,^;0) + 0(r,2 



7=7o+A7 



(56) 



where A7 is linear in ro, and we used the fact that 4(0, 0) = 0. 

If we square the amplitude in Eq. (56), we can calculate ai at 7 near 70 as 

2 



5127r^ 



;fto(0,^;0) 
d7 V3 



7=70 



A7 + rofi(0,^;0) 



(57) 



where we used the fact that = at 7 = 70. The leading term here is quadratic in tq 
and purely determined by up-to-NLO scattering amplitudes. Higher-order corrections to 
the scattering amplitude only start to affect a at order Tq. In other words, even though the 
NLO a near 70 is in the order of Tq, it is purely determined by LO and NLO t-matrices. 



From Eq. (25) and Eq. (28), we derive 



d ~. 27 , 



:^-/f„(0,^;0) 



7=70 



7=70 



and 



f.(0,^;0) 

where we again applied i^o(0, 0) = 0. 
Therefore, ai near 70 is 



^1(0, ^;0) 



1 _ ,-15lo (tio 270. 



0) 



(58) 



(59) 



5127r^ 



1^^0(0, ^;0) 



27 . 



7=70 



A7 + roiri(0,^;0) 



V3 



m 



i-z4^iro(^,^;0) 



(60) 



The next-to-leading-order recombination minimum ai = is thus determined by 



d 



d7 



27 



— /To 0,-L;0 



7=70 



A7o + roiei (0,^;0) = 0, 



(61) 



which leads to an NLO shift in the position of the recombination minimum of: 

^1(0, ^;0) 



A70 = -ro 



£^o(0,^;0) 



(62) 



7=70 
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D. Atom-Dimer Resonance 



The atom-dimer scattering length diverges for positive scattering lengths for which B = 
—'-f^/m. At this value of the scattering length, a three-body state lies exactly at the atom- 
diatom threshold. This feature shows up as resonant behavior in the atom-dimer relaxation 
rate, a process in which shallow dimers are transfered through collision processes into deep 
dimers. The two-body binding momenta for which these resonances occur are denoted (at 
LO) by 7*. They thus obey the relation 

BoM = --. (63) 

m 

Discrete scale invariance in the leading-order bound-state spectrum implies that if 7=,, is a 
solution of Eq. (63), then so are the quantities e'^'^^^^j^. The scale invariance is softly broken 
by To /a corrections. Here we will calculate the NLO corrections to these 7*'s. 

We will assume that the position of a particular resonance is shifted to 7=,, + A7* at NLO. 
The three-body binding energy at the atom-dimer threshold must then obey up to terms of 
relative order C(7^rQ): 

5o(7* + A7,) + 5i(7, + A7,) = J^* + ^^*)\ (64) 

m 



We now expand both sides of (64) in powers of ro and retain only terms up to O{ro): 

5o(7,) + — 7=7,A7, + 5i(7,) = 2 . (65) 

(27 mm 



Using Eq. (63) we find a NLO correction to 7* 



mfii(7*) 

27* +"^^?Pl7=7. 



We note that Eq. (43) implies that A7=i, is linear in tq. 



However, calculating A7^, according to Eq. (66) results in numerical difficulties, because 
-81(7*) and the denominator are both zero to within the numerical accuracy of our calcula- 
tion. 

The fact that the denominator should go to zero is clear from the expression given for 
Boil) in [13]: 

K = —if sin ,^ 
7 = HcosC, 

H = (e--/^«)"-"*«,exp[A(0/(2so)], (67) 
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where k = y/—mBo. When 7 is near 7=,, the function A(,^) can be expanded in power of 
(— 7r/4 — ^)^/^. For completeness we give the expression of Ref. [T3|, with the coefficients 
that were determined numerically there: 



A = 6.04 - 9.63(-^ - 0'/' + 3.10(-^ - 0, 



(68) 



where ^ = — 7r/4 corresponds to the three-body bound state crossing the atom-dimer thresh- 



old at 7 = 7^,. We calculate dn/d'y from Eq. (67) and find that at ^ = — 7r/4 



dK, 
d'y 



7=7* 



dH , TT 



dA. 
d^ 



2so 



^=-,r/4 d( ^ 



(69) 



?=-7r/4 



since — )■ cxd when ^ — )■ —n/A. dBo/d^ at the atom-dimer threshold is therefore 



dBn 



m- 



d'j 



-2k 



7=7* 



dn 
d'j 



-27.. 



(70) 



7=7* 



This indicates that the denominator in Eq. (66) indeed goes to zero at the point of interest, 



which causes an accuracy problem in numerically calculating A7* from the bound-state side. 

We calculate therefore instead A7=i, from the scattering amplitude, i.e. evaluate the atom- 
dimer K-matrix as a function of 7 along the threshold line E = —7^. The LO three-body 
scattering length af^ is related to Kq at this energy by 



4°^ = -?i^o(0,0;-f). 



(71) 



Now af^ — )■ 00 at 7 = 7*, and so the on-shell Kq has a pole of order one at 7 = 7*: 



7-7* 



+ 7^o(7)• 



(72) 



Similarly the NLO shift of the position of the atom-dimer resonance A7=i,, can be written as 

^ad ^ ^ad 



Ko{0, 0; -7') + 7^oi^o(0, 0; -7') + roK,{0, 0; -7' 



ryad 
^0 



+ 



7-7* 7-7* 



A7, 



7 - 7* - A7* 



(73) 



The shift A7=i, is therefore calculated as: 

lim^_^,(7-7,)2i^^i(0,0;-72) 



A7* = ro- 



(74) 
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Our numerical calculation shows that both the numerator and the denominator in this form 
are finite, and so calculating A7^, by evaluating the atom-dimer K-matrix at different 7's 
along the threshold line is an accurate procedure. 



However, we still have to show that Eqs. (66) and (74) are equivalent for calculating 



A7*. The bound-state form of the leading-order atom-dimer amplitude near a bound state 

mZo(0,0) 



of energy -80(7) (Eq- (36)) is: 
to(0,0;-7') 



-7^-50(7)+''° 
mZo(0,0) 



(7 - 7*) 



27* + 



7=7* 



(regular). 



(75) 



This relates Zq in Eq. (1361) to Zg"'^ in Eq. (1721) as 



Zo(0,0) 



27* + 



dBn 



7=7* 



(76) 



7=7. 



Similarly we expand (|38|) about 7^, and so relate ZqBi to the numerator in Eq. (|74|) as: 

dBo 



ZqBi 



27* + 



7=7* 



6/7 



7=7*, 



ro lim (7-7,)^K,(0,0;-7^), 

7->7* 



(77) 



which explains why i?i — )► as 7 — )■ 7=,, , and, moreover, shows that the coefficient of this zero 
is precisely what is needed to render the expressions obtained for A7^, from the bound-state 
and scattering-state side equivalent. 



E. Three-atom resonance 



When a state in the three-body bound-state spectrum crosses the zero-energy threshold 
at negative scattering length, three free atoms can form a zero-energy trimer state. This 
phenomena is called a three- atom resonance, and results in a maximum in the three- atom 
recombination rate. It occurs at a value of 7 denoted by 7_. (Since 7_ < this does not 
correspond to the binding momentum of a dimer, but it is still the inverse of the atom-atom 
scattering length where this feature occurs.) At leading order the condition for this three- 
atom resonance is i?o(7-) = 0. Furthermore, discrete scale invariance of the leading-order 
bound-state spectrum then results in values of 7_ being related by the universal scaling factor 
g7r/so_ (Here and below 7_ denotes the leading-order position of the three-atom resonance.) 
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The NLO correction to 7_ is found at the NLO zero-energy threshold {Bq + i?i)(7_ + 
A7_) = 0. We thus have: 

dBoi^) 



Boil-) + 



A7_ +5i(7_) 



0. 



7=7- 



As i?o(7-) = is given at LO, we find the NLO correction to 7_ 

Bid-) 



A7_ 



dBoi-y) 






7=7- 



(78) 



(79) 



and so A7_ is hnear in tq. Our numerical studies show that neither the numerator nor the 
denominator in this equation are equal to zero at 7 = 7_. Our result is: 

dBoi-f) 



-0.9847_, 



(80) 



7=7- 



where the coefficient is calculated with a numerical accuracy of about 10 ^. This is in 



contradiction to the form provided in Ref. [13], which predicts 
Appendix IbI) 



(7) 



-27_. (See 



7=7- 



V. THE SUBLEADING THREE-BODY FORCE AT NLO 



We will show in this section explicitly that the NLO counterterm contains a scattering- 
length-dependent piece that will require a second experimental datum for renormalization 



if scattering-length-dependent processes are considered. To do this we reconsider Eq. (43), 
the expression for the NLO shift to the binding energy: 



Bi 



IT 



d<i <f^^^ '"!'^~"% n<D 



-7 + V3gV4 - tuBq 



8gi(A)ro 
(7rA)2 



n 2 



:r(g) 



(81) 



-7 + \/3gV4 - niBo 

We will use the divergence structure of this observable to determine the behavior of Hi{A) as 
a function of A, up to corrections ~ 1/A. The divergence structure of any other observable 
computed to NLO will be similar, and so it suffices to perform this calculation for Bi. In 
particular, we will expand both the explicit integrals and the behavior of Hi{A), in powers 
of A, and demand that the linear-in-A and log(A) divergences cancel. 

In order to perform this analysis we need to know the large-momentum behavior of each 



term in Eq. (81 ). At large momenta q the function r(g) is known in the form of an expansion 

20 



in powers of 7/g (see Appendix of Ref. [6] where we have corrected an error in the result for 

(82) 



r(g) oc — + — + .. . , 
q 



with 



Zq 

Zl 



sin ^so In -2 
2 

7! 



where 



and 



C_i 



lis) 



C_i I sin ^so ^ + ^^S ^-1) 
/(iso - 1) 



l-/(zso-l) 



I sm 



Vs. 



s cos 



(83) 
(84) 

(85) 
(86) 



Inserting the asymptotic form of r(g), (82), up to ~ l/?^, into Eq. (81), and evaluating 
the first integral shows that Hi has to absorb both a linear divergence and a logarithmic 
divergence proportional to 7. In order to cancel these cutoff dependencies we will thus write 
Hi as 

ifi = A/iio(A)+7/^ii(A). (87) 
Analytic expressions for both /iio(A) and /in (A) can then be obtained by inserting the 



expansions ml\ and (82) in Eq. (81), while also expanding all explicit functions of q in 



powers of 7/g. In this way we find: 



rfg gM 1 + 



47 



2^27 

Zq H ZqZi 

q 



8A 



7r2A2 
1 

TT 
8 

'^A 



{hw + '-j^hii 



2 /■ , / 27 

Til '"l' + vr,y, 



^0 + -Zl 



2 47 2 27 

Zq + ~^^0 ^OZl 

V3g q 
(hio + j^hiij ^ j dq (^zq + 



27 , 7 

—7=-Zo + -2i 

V3g q 



+ . . 



(88) 



where ( is finite, the dots represent finite parts of the integration, and defines an integral 
that is regulated in the ultraviolet by a cutoff A and whose infrared regularization (if any) 
is unspecified. 
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In order to simplify the notation we denote integrals that contain a product of the z 
functions by a W. The first and second indices indicate the 2;-functions in the integrand, 
while the third index gives the power of q that resides in the denominator of the integrand, 
so 

m^n^-tdq'-^. (89) 

In the same spirit of notational convenience and compactness we define: 

1 /-^ . 



Zrnn=- dq ^ . (90) 

All integrals W/mn and Z^n can be evaluated analytically, and this is done in Appendix |Aj 



The divergences linear in A in Eq. (88) are then cancelled by requiring: 



>Vooo(A) + ^^oo(A) = . (91) 
Meanwhile, the divergence which is logarithmic in the cutoff is canceled by 
-|wooi(A) + Won(A) + ^^o'o(A) 

+ ^Zoo(A) (^^^01 (A) + 2ii(A)j = . (92) 



By using the results for these integrals which are given in Appendix |A[ we derive the 
analytic forms for the NLO piece of the three-body force. First, 



, ... 37r(l + s2) v^l + Asl - cos (2so ln(A/A) - arctan 2so) , , 

64^1 + 4sg sin^ (so ln(A/A) - arctan So) 

Since A is determined by the LO renormalization condition (e.g. A = 13. Ik,,, in the unitary 



limit), Eq. (93) is a prediction for hiQ^K). Meanwhile, Eq. (92) gives 



y37r(l + sg) (l + |C_i|cos(argC_i)) , 

/iii(A) = ^ (M-K\ 1 s^'^K^IN 

io sm (so ln(A/Aj — arctan So j 

v^7r(l + si) sin (2so ln(A/A)) + |C_i | sin (2so ln(A/A) + arg C_i) 



32so sin^ (soln(A/A) — arctan so) 

v^7r(l + s2)3/2 cos (soln(A/A)) + cos (soln(A/A) + argC_i) 



16so sin^ (so ln(A/A) — arctan sq) 

cos (2soln(A/A) — arctan(2so)) 

Vl + 4s2 



X 



(94) 



^ The result for, e.g. Zqq{K), would seem to neglect the effect of the infrared regularization, which could 
affect the answer for hn. However, the combination of the infrared-regularization dependence of Zqi[K) 
and Zii{K) is zero, and those of >Vooi(A) and Won (A) can be absorbed into the numerically fitted 
parameter /Lt, which is included in /in (A). 



22 



The /i in Eq. (94) subsumes information on the finite part of (9(A''), and its value is de- 
termined by the renormahzation conditions at NLO. Numerically this piece is of the same 
order as In A, as long as A is not extremely large. 

We compare the analytical and numerical results for a particular renormahzation condi- 
tion in Figs. |5] and [6j Note that in order to do this comparison in Fig. [6] we have fitted the 



value of to the numerical results. However, the coefficient of the log(A/yu) term in Eq. (94) 
is still predictive, and is confirmed by comparison to the numerical results. The analytic 
and numerical results agree to within our numerical accuracy of 10^^ when /i = 0.990k^,. 

The renormahzation condition chosen here is k^, = 1 at LO and NLO, and 70 = k^,/0.316 
at NLO. This corresponds to a recombination minimum that receives no range corrections at 
NLO. The value of /i is affected by the choice of NLO renormahzation conditions. We have 
checked that the numerical results for hu (A) obtained with other renormahzation conditions 
can be well described by choosing alternative values of /i. 

In both figures we plot the inverse l//iio and l/hn of the low-energy coefficients, and so 



the bumps that can be seen in Fig. p^are the result of zeros in Eq. (94). The competition 



between the different terms in Eq. ([94j) is clearly seen in the comparison. In particular, the 
double pole has a coefficient which grows with In A, and this produces a turning point which 
is in increasing proximity to the position of the dominant triple pole as A grows. 

In principle, the details of the regularization at NLO can affect the discrepancy between 
analytical expressions and numerical results in the NLO three-body force. In order to 
eliminate this effect, we first calculated the LO amplitude, Iq, at a large cutoff, A = IO^^k^,. 
Then we insert the resulted to into NLO calculations integrated to a relatively smaller cutoff, 
A < lO^K*, which suppresses the corrections from the details of the regularization to < 10^^. 



VI. UNIVERSAL RELATIONS AT NLO 

In the limit ro/a = 0, the two-body scattering lengths at which key experimental features 
occur, such as the atom-dimer resonance at positive scattering length a = a^,, the three-atom 
resonance at negative scattering length a = a_, and the recombination minimum at a = oq, 
obey a set of universal relations. Indeed, all three features can be related to the binding 
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10^ 10^ 10^ lO"^ 10^ 10^ 



FIG. 5: /iio(A): The dots are numerical results for hio with the NLO calculation renormalized such 
that the LO prediction of k=k = 1 (A = 13.1) is maintained. The solid line (red) is the analytic 
function /iio(A) given by Eq. (|93|) with the same parameter A. 



momentum of the three-atom state in the unitary limit, [13 



a* = 0.0708^;^ , (95) 
a, = -1.51k;' , (96) 
ao = 0.316k;^ . (97) 

We can eliminate in the equations above and then rewrite them in terms of the inverse 
scattering length (7 = 1/a at LO) 

7* = O^-fo , 7- = O-lo , 1^* = 6*0070 , (98) 

where 9^ = 4.47, 9^ = -0.210 and 9^ = 0.316. 

Once the two-body effective range, tq, is non-zero we can extend the LO universal relations 



such as Eq. (98) to include its effects. However, as shown above, these relations will now 



require two three-body parameters as input to predict a third one. This means that, before 
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10^ 10^ 10^ lO"^ 10^ 10^ 

FIG. 6: /iii(A): The dots are results of a numerical calculation of hu with the additional NLO 
counterterm fitted to maintain the value of 70 = k*/0.316 predicted at LO. The LO renormalization 
condition, = 1 (A = 13.1), is also kept at NLO. The solid line (red) is the analytic function 
/ill (A) given by Eq. (|94]), with fitted to 0.99k* in order to obtain good agreement with the 
numerical values. 

predictions can be made, two three-body observables have to be determined, either from 
experiment or from the underlying microscopic interaction that governs two- and three- 
body properties. To illustrate this, we assume that the positions of the two features 7=,, and 
7o are known. The effects of the effective range perturbation can then be included by adding 
a term linear in ro to the LO universal relation between 7=,, and 70 

7* = 6^*70 + ^"0^^*70 ■ (99) 

The real and dimensionless number Z,, is non-universal in the sense that it depends on the 
renormalization condition, and thus on the type of physical system under consideration. 
The universal relations between other features are modified in the same way. For example 

A7_=roX_7o, (100) 
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where 7_ can be predicted if the relation between Z,, and X_ is known. Similarly, the shift 
from the LO value of be written as 



Ak* = roXoo7o • 



(101) 



Note that in Eqs. (99)-(101) we have employed 70 to make up the dimensions of the NLO 



correction. This is a matter of convenience, since all dimensional quantities are related at 



LO by formulae such as Eqs. (95)-(97). We also used the fact that the NLO shifts in these 



three-body observables are strictly linear in tq in our perturbative analysis. 

While the parameters that we have introduced above are not universal, the relations 
between them are. These are linear as a consequence of first-order perturbation theory. We 
have obtained these linear relations numerically and verified their stability for a wide range 
of possible renormalization conditions. In the case of the renormalization scheme discussed 
in the previous paragraph (70 fixed at LO and NLO, and 7* the additional input at NLO 
needed to fix hu) we found 



(102) 



with 



^00,* 

V-,* 



0.256 , 
-0.286 , 
1.65 X 10"^ 
-2.04 X 10" 



(103) 



Combining Eqs. (99), (100), (101) and (102) with each other leads to relations that contain 



only observable features and the parameters in Eq. ( 103 ) 



7- = ^-7o + ^-,*?^o(a_)7o +^-,*^^7^ (7* - ^*7o) 



2 , '"o(aoo) 

00,* / \ 

ro(a*) 



(7* - 6**70) 



(104a) 
(104b) 



It is important to note at this point, that different universal relations can be generated 
depending on which three-body features are assumed to be known. 



26 



We can now apply our results to recent experimental measurements of recombination in 
ultracold gases. Litliium-7 atoms in the \F = 1 mp = 0) hyperfine state are one system 
for which the scattering length and effective range are known as functions of the magnetic 
field. Gross et al. [23] measured a number of recombination features for this system. On the 
positive scattering length side of the resonance, they were able to measure a recombination 
minimum at oq ~ HGOob and the position of an atom-dimer resonance at a^, ^ 2900^. 
On the negative scattering length side, they measured a recombination maximum a[_ ^ ^ 
—2640^. The superscript indicates that this a_ represents the point where the next 
deeper trimer branch relative to Oq and a^, crosses the three-atom threshold. Ref. [23] also 
contains the results of a calculation for the magnetic field dependence of scattering length 
a and effective range tq. We will employ these results to illustrate the use of the relations 
presented above. 

The values of the effective range at the scattering lengths where experimental features 
occur are 

roiao) = -M.baB 
ro{a*) = -7A.7aB 
ro (aL^^) = 27.2aB, 

We now use the values of two of these observables to predict the third one (whose measure- 
ment was already reported in Ref. [23])- In particular, we will predict a^, after renormalizing 
to and a^~\ In doing this we assume that the low-energy constants Ho, hiQ, and hu de- 
pend only weakly on the magnetic field and are the same for all three experimental features 
we consider. We start out again by relating the known features through the LO universal 
relations and a shift linear in the effective range times a non-universal number that contains 
the information about the NLO renormalization condition 

70 = -0.2107!"^ + ro(ao) ■ X^-^i")' , (105a) 
7, = -0.9397i~^ +ro(a,) ■Xi-)7i"^' , (105b) 

We again can numerically extract the universal parameters that relate the parameters Xq""* 
and xi""* via 

X(-) = -0.309 + 7.17Xi"^ , (106) 
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here we haven chosen to write out the numerical values of the universal parameters directly 
to avoid notational clutter. We can then combine these equations into one universal relation 
that predicts the position of a^. 

The corresponding universal relation is given by 

^2 



% = eth^^ - 0.309 ro(a.)7i-)' + 7.17^ (70 - ^^7^^ 



(107) 



where the parameters of the universal LO relations are given by = -0.939 and O'q 
-0.210. 

We can alternatively generate a relation that gives the NLO result for 7^, as the relative 
shift to the LO universal equation between 7^, and 70. 



7* 



- 7.02 ro(a,)7^ + 0.566 



rpja*) / (-) 



7- 



108) 



where 9, = 4.47, = l/^J"^ 



ro(a^^) 

4.75 and we have again immediately inserted the numer 



ical values of the NLO universal parameters. Equations (107) and (108) will generally give 



different numerical results for 7* which provides a lower bound on higher order corrections. 

Using the numerical values for the positions of three-body features and the numerical 
values for the effective ranges, we obtain 



a, = (271 - 108 + ...)aB 
a, = (257- 2 + ...)aB , 



(109) 
(110) 



from Eq. (107) and Eq. (108), respectively. 



Our predictions of a^, in both renormalization schemes are consistent with our previous 
findings [Pj. Combining both results and including a conservative estimate of higher order 
corrections, we obtain therefore the same result as quoted in Ref. [9] 



a, = (210 ± U)aB . 



;iii) 



VII. CONCLUSION 



In this work we have presented a perturbative calculation of next-to-leading order (in 
£/|a|) corrections to universal three-body physics. We have shown that an additional three- 
body counterterm is required for renormalization if and only if scattering-length-dependent 
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quantities are considered. The inverse scattering length therefore plays a similar role in the 
SREFT to that of the pion mass in chiral perturbation theory. Counterterms proportional 
to the inverse of the scattering length occur at orders beyond leading, and these must be 
fitted by considering scattering-length-dependent data. 

The advantages of the perturbative analysis we have presented here are twofold. First, 
it allowed us to derive analytic expressions for the next-to-leading-order shifts of resonance 
positions in recombination experiments. Second, it permitted an explicit treatment of the 
renormalization of divergent integrals while keeping the LO counterterm fixed. This pro- 
duced analytic forms for the running of the NLO pieces of the three-body force. 

Our analysis is applicable to systems for which the scattering length and effective range 
are known as a function of the magnetic field. If data on three-body processes at different 
values of the two-body scattering length exists then effective-range corrections to recombi- 
nation features can be treated in the manner described above. Data from the Bar-Ilan and 
Rice groups [231 121! Lithium-7 recombination were recently analyzed in this way [0]. 

A calculation of the divergence structure of bosonic observables in SREFT at next-to- 
next-to-leading order C(£^/a^) is underway [25] . 
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Appendix A: Relevant Integrals 

The following functions are used in the derivation of the analytical expression for the 
NLO three-body counterterms hiQ and hu. 

1/"^, ^ , . f . <1\ A ./A 

^00 = — = — ag sm So In ^ = ^^^= sm Sq In ^ — arctan Sq , (All 

Try TT J V A/ n^yl + sl V A / 
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txsq \ A 
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: — / — sm So m ^ + argC_i 

71 J q V A 
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-^=— cos 1^ So In ^ + arg C_i 
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Integrals that contain a product of the z functions are denoted by W. The first and 
second index indicates the order of the two 2;-functions in the integrand's numerator, and 
the third index gives power of q in its denominator. I.e. 
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Appendix B: Parameterization of the three-body spectrum on the negative- 
scattering-length side 



We foUow Ref. [I3] and define Bq on the negative 7 side as 



5o(7) = -l^Q 



1_ 

7- 



(Bl) 



where the dimensionless function Q{1) = at threshold. By substituting this into Eq. (79) 
we have 



A7_ 



Bid-) 



(B2) 
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The calculation of A7_ is numerically accurate because G'{1) 7^ 0. We find G'{1) = 0.98. 
That 7^ is supported by the empirical fitting formula in [13] near 7 = 7., where k 



and 7 obey Eq. (67), while the function A(^) obeys a different expression from Eq. (68), 
because now we are near the three- atom resonance, not the atom-dimer resonance, which is 
indicated by ^ = — vr: 



dBo/d'y in this region is 



-7r, — 



: A = -0.89 + 0.282 + 0.25;z^ 
^ = (vr + e)'exp[-l/(7r + e)1 



(B3) 



^5^_2i7sine(2i^f sine + cose) 



d'-f 



if cose -sine 



(B4) 



whose numerator and denominator both go to zero when C, — > — vr. After applying I'Hopital's 
rule we calculate dBo/d'j at 7_: 



dBo 



d'j 



2H 



7=7- 
32 A /^t-a 



1 



1 d^A 
2so (1(2 



(B5) 



Equation (B3) indicates that d A/d^ = at ^ = — tt, therefore 

dBo 



d^y 



2H = -27_ 



(B6) 



7=7- 



and so 



G'il) = 2. (B7) 
Thus, the function G'{1) is finite and nonzero from the empirical fitting formula in [T3] . 



However, the value in Eq. I\B7« is larger than the correct numerical value in Eq. (B7) 
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